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ABSTRACT 



Alms. We suggest an explanation for the high-frequency power excess surrounding active regions known as seismic halos. 
Methods. We use numerical simulations of magneto-acoustic wave propagation in magnetostatic sunspot model. 
Results. We propose that seismic halos can be caused by the additional energy injected by high-frequency fast mode 
waves refracted in the higher atmosphere due to the rapid increase of the Alfven speed. Our model qualitatively explains 
the magnitude of the halo and allows to make some predictions of its behavior that can be checked in future observations. 
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1. Introduction 

Almost since the discovery of the 5-min solar oscillations it 
is well known that the oscillation power is re duced by some 



40-60% in the p hotos pheres of suns pot s (Lit es et al. 
'Abdelatif et al.', '1986; Brown et alj, Il992i : iTitle et al 



1982; 



19921 : 



.Hindman fc Brown, 1998). Later it was found that the high- 
frequency non-trapped wave power shows a suspicious en- 
hancement in rings surrounding active regions, both in the 
photosphere (iBrown et al.L 1992h and in the chromosphere 
(jBraun et al.l . Il992l : iToner fc Labontd 11993). These power 
enhancements are known as "halos". Their observational 
properties can be summarized as: 

(i) The power enhancement is observed at high frequen- 
cies between 5.5 and 7.5 mHz for waves that are usually 
non-trapped in the non-magnetic quiet Sun. 

(ii) The acoustic power measured in halos is 
larger than in the nearby qu ite Sun by abo u t 40- 
60% ([H indman & Brown, 1998; lBrau_n fc Lindsevl. Il999 : 
Donea et al.. . .200ft : iJain fc Haberl . 20021 : iNagashima et al 



2007D . 



(iii) The halos are observed at intermediate longitudi- 
nal magnetic fluxes {B) = 50 — 300 G, while the acous- 
tic po wer is usually reduc e d at a ll frequencies at larger 
fluxe s f Hindman fc Brown], Il998l : iThomas fc Stanchfleldl . 
120001 : [jain fc Ha ber. 20M)- 

(iv) The radius of the halo increases with height. In 
the photosphere the halos are located at the edges of ac- 
tive regions, while in the chromosphe re they extend to a 
large portion of the nearby quiet Sun (iBrown et al.L 1 19921 : 
iBraun et all 1 19921 : Pfh omas fc Stanchfleld l. 12000^. 

(v) The power increase in the halo is qualitatively sim- 
ilar in sunspots, pores and plages. 

(vi) Signiflcant reflection of the upcoming acoustic ra- 
diation at 5-6 mHz is detected in active regions, unlike 
the behavior of such high -frequency waves in the quiet Sun 
(jBraun fc Lindsevl . 12000 ) . 



While several plausible mechanisms have been pro- 
posed to explain the acoustic power reduction for the 
stronges t flelds in active regio ns (e.g. MHD mode con- 
version; I C ally fc BogdanI , Il997f ). no accepted theory ex- 
ists to explain the power enhancement in acoustic halos. 
The increase of the high-frequen c y aco ustic emission, ini- 
tially proposed by i Brown et al.l (|l992l ) and iBraun et al.l 
(|l992f ). seems not to flnd observational conflrmations, 
since the observed continuum intensity does not show the 
halo effect (Hindman fc Brown, 1998; Jain fc Haber, 20Q2|). 
Alternatively, the latter authors propose that the velocity in 
the surroundings of active regions may become fleld-aligned 
and some type of incompressible waves may be responsible 
for halos. This, however , lacks any observational evidence. 
Recently. i Kuridze et al.l ^2008) suggested yet another mech- 
anism based on acoustic waves trapped in fleld-free atmo- 
spheres lying below small-scale magnetic canopies of net- 
work cores and active regions. Interestingly, halos were ob- 
served recentl y in MHD simulations of waves in magnetic 
structures bv iHanasogd C2 008) and S hely^ g et al. (200l). 
Based on his simulations, iHanasogd (|2009l ) suggests that 
the power enhancement in halos is due to magnetic fleld in- 
duced mode mixing resulting in preferential scattering from 
low to high wave numbers. 

In this Letter we propose a mechanism based on the 
fast MHD mode refraction in the vicinity of the transfor- 
mation layer (where the Alfven speed va is equal to the 
sound speed cs) that is capable to explain several obser- 
vational properties of halos. In addition, we predict some 
new properties that can be obtained from observations in 
the future to conflrm or discard this explanation. 



2. Description of methods 

We perform 2D numer ical e xperiments that are essen- 
tially similar to those by iKho menko et al. (20091), to study 
the adiabatic propagation of magneto-acoustic waves ex- 
cited by a single source located at sub-photospheric lay- 
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Fig. 1. Snapshots of the vertical velocity (scaled with the factor y/p) at elapsed time 60 minutes in the simulations with 
harmonic source placed in the sunspot model (left panels) and in the quiet Sun model (right panels). Panels from top to 
bottom: source periods 100, 200 and 300 s. The upper plots on each panel give the r.m.s. amplitude distribution at the 
line formation level. The curve from the sunspot simulations is repeated in dashed line on the right panels. The contours 
of constant c|/'U^ are marked with numbers. The dashed-line contours indicate the "line formation" heights in both 
models. Dotted lines: height where the period of the source coincides with the cut-off period. The black inclined lines are 
magnetic field lines. 



ers of a magneto-static sunspot model. The numerical 
MHD c ode is described in deta il in iKhomenko fc ColladosI 
(|2006[ ): IKhomenko et all (|2008[ ). T he unperturbed magne- 
tostatic sunspot model is taken from'Khomenko fc ColladosI 
(jloOS). The simulation domain has 80 x 10.5 Mm^ in hor- 
izontal and vertical directions, respectively, with a resolu- 
tion of dx = 0.15 Mm and dz = 0.05 Mm. The whole do- 
main contains magnetic field, but it becomes weak and 
dynamically unimportant further than ~20 Mm from the 
sunspot axis. The maximum field strength in the photo- 
sphere is around 1 kG. Our sunspot model has a Wilson 
depressio n. With the help of the SIR radiative trans- 
fer code (|Ruiz Cobo fc del Toro Iniestal , Il992l ), from the 
known distribution of thermodynamic parameters in ge- 
ometrical height, we calculated the optical depth scale, 
logTs, for each horizontal point of the MHS sunspot model. 
The photospheric level defined by the optical depth scale 
logTs = is located 300 km deeper at the sunspot axis 
compared t o its location 40 Mm far from the axis (see Figs. 
1 and 2 in IKhomenko et al.L 120091 ) . We define two refer- 
ence levels of optical depth: logrs = ( "photosphere" ) and 
logTs = —1.6 ("line formation"). At a horizontal distance 
X = 40 Mm from the sunspot axis, the "photosphere" is 
located 500 km below the top boundary of our simulation 
domain. The "line formation" level is located 160 km above 
the "photosphere". At 20 Mm from the axis, the va = cs 
level is above the "line formation" level by about 300 km, 
and at the axis it is located some 200 km below. 

In the first set of experiments the source is placed at 
three different horizontal distances Xq = 20, 30 and 35 
Mm from the sunspot axis and at Zq = —700 km below 



the photosphere. The tempora l behavior of the source is 
descr i bed by a Kicker wavele t ([Parchevskv fc Kosovichevl 
I2OO9I : IKhomenko et al.L l2009l ) with a central frequency of 
3.3 mHz. In the second set of experiments the source is 
placed at Xq = 30 Mm and Zq = —450 km, and it is har- 
monic and continuous in time with three different frequen- 
cies = 10, 5 and 3.3 mHz. In all the cases the sunspot 
simulations are accompanied by non-magnetic simulations 
in the modified mo del S of 'Christensen-Dalsgaard e t al.l 
(1996) (see Parchevs ky fc Ko sovichev, 2007) with exactly 
the same properties of the source and numerical treatment. 
The duration of simulations is about two physical hours. 



3. Results 

Fig. [1] shows snapshots of the vertical velocity in simula- 
tions with the harmonic source. When the wave frequency 
is well above the cut-off frequency {u = 10 mHz) the waves 
are propagating in the quiet Sun (Fig.[T)3). This situation is 
different in the sunspot model. Significant refiections can be 
appreciated in Fig. [1^ visible as an interference wave pat- 
tern around X = — 20 : — 10 Mm. The reflection of the high- 
frequency waves is produced in the vicinity of the trans- 
formation Va = Cs layer. The fast (acoustic-like) waves 
generated by the source are transmitted as fast (magnetic) 
waves in the upper atmosphere where the Alfven speed 
is larger that the sound speed {va > cs). In the mag- 
netically do minated layers these waves are refracted (see 
Khomenko fc Collado s. 2006) and are returned back to the 
sub-photospheric layers where they interfere with the waves 
coming from the source. This behavior of waves is clearly 
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Fig. 3. Ratio between the photospheric wave power in the sunspot model relative to the quiet Sun model ((Pspot — 
Pquiet)/Pquiet) cis a fuuctiou of horizontal position and frequency for three simulations with the wavelet source located at 
Xq = —35 (a), —30 (b) and —20 (c) Mm. The color scale at the bottom panels ranges from -50 to 50% (red means power 
excess). The upper panels show the power ratio averaged over the low-frequency and high-frequency bands (marked on 
the figure). The vertical dashed line located at X = marks the sunspot axis. Another dashed line marks the location 
where va = Cg in the photosphere. 
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Fig. 2. Wave paths of the fast mode launched from the 
lower turning point at X = —30 Mm, Z = —1 Mm prop- 
agating through the sunspot model for frequencies 7.5, 5.5 
and 3.5 mHz (from top to bottom). For clearness only the 
upper part of the model is shown (not to scale). Blue in- 
clined lines are magnetic field lines. The yellow contours 
mark the layer where the wave frequency is equal to the 
cut-off frequency. The white contours are va = cs- The 
background image is logio('^A)- Each white dot on the tra- 
jectory is separated 3 minutes in time. The red dashed lines 
are the "photosphere" and the "line formation" levels. 

seen on the movie of this simulation, attached as on-line 
material to this paper. The r.m.s. vertical velocity ampli- 
tude measured at the line formation level (dashed curves 
in Fig. [1^ and b) shows a suspicious bump on the left from 
the location where the logrs = —1.6 contour crosses the 



= contour (see upper plots on each panel). This bump 
is absent in the quiet Sun simulation. We propose that the 
increase in the high-frequency power in halos surrounding 
active regions can be produced by the additional energy in- 
jected by the fast mode waves refracted in the magnetically 
dominated layers back to sub-photospheric layers. The pres- 
ence of the upward and downward propagating wave energy 
manifest itself as the wave interference pattern. Note that 
this mechanism does not necessarily imply that the high- 
frequency fast waves are trapped in sunspot s, since a part 
of their energy can leak into the slow mode waves after 
each mode transformation at the va = cs layer. However, 
this mechanism produces significant reflections of the high- 
frequency waves that otherwise are propagating in the quiet 
Sun. 

When the wave frequency is below the photospheric cut- 
off frequency {v = 3.3 mHz, Fig.[T^, f) the fast mode waves 
become evanescent before reaching the transformation va = 
Cs layer. In this case neither refraction nor interference can 
be produced because of the absence of fast to fast mode 
transmission. Fig. [1^, f (upper plots) shows that the r.m.s. 
amplitude distributions are very similar in the sunspot and 
the quiet Sun simulations and that the amplitudes in the 
sunspot case are always lower than in the quiet Sun. 

The u = 5 mHz case (Fig. [Th, d) shows an intermediate 
situation, where the waves are reflected both due to cut- 
off effects and the magnetic effects. In this case the cut-off 
height (dotted lines) almost coincides with the va = cs 
height. Some power excess is still present in the sunspot 
case compared to the quiet Sun. The bump is weaker and is 
now located more near the source as the interference hap- 
pens at another location due to the different wavelength. 
Note that the r.m.s. velocity distribution is less smooth for 
the V = b and 3.3 mHz compared to = 10 mHz, as several 
nodes are present in the horizontal (and vertical) direction 
due to the evanescent character of waves. 

The eikonal solution for the fast mode wave 
120061 : iMoradi Callvl . [20081 : iKhomenko et al.l . [2009h allows 
to get a more complete physical picture of the wave be- 
havior at different frequencies. Fig. Ogives the wave paths 
of the fast mode waves launched from their lower turning 
point (that roughly coincides with the location of the source 
in the simulations). At high frequencies (5.5 — 7.5 mHz) the 
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waves penetrate higher up in the atmosphere above the hne 
formation layer where they are refracted back down due 
to the presence of the magnetic field. In the low-frequency 
case (3.5 mHz, Fig. [St) the waves are sharply reflected from 
the cut-off layer and are less affected by the magnetic field. 
Note that the power excess (halos) at high-frequency should 
form at distances where the refraction of the fast mode oc- 
curs above the line formation layer, i.e. where the va = cs 
layer lies above the photosphere. Otherwise it would not be 
detected in spectral observations. This happens in regions 
of intermediate field strengths. Under no circumstances can 
the halo form in the "umbral" zone with strong fields where 
Va is larger than cs in the photosphere. 

A more realistic situation with the source emitting a 
spectrum of waves is considered in Fig. [3l In this figure 
we compare the wavelet source simulations with different 
source locations. We represent the ratio between the pho- 
tospheric Fourier power in the sunspot and in the quiet 
Sun at different frequencies as a function of horizontal dis- 
tance. We conclude that, independently of the position of 
the source, the overall picture is very similar especially for 
the sources at —35 and —30 Mm. For the source located 
at the —20 Mm from the axis, the magnetic field plays al- 
ready some role and modifies the properties of the source 
compared to the purely non- magnetic case, thus contam- 
inating the detailed picture. Still, even in this case, halo 
power increase is present around —10 Mm, being weaker 
than in the other two cases. Fig. [3] shows a power excess at 
high frequencies above 5.5 mHz at distances —20 : —10 Mm 
from the sunspot axis (the latter position coincides with 
the location where va = cs in the photosphere). At low 
frequencies, a power deficit is observed in the sunspot. The 
power deficit is also present at high frequencies in the "um- 
bral" region. The power averaged over both, low-frequency 
and high-frequency, bands shows a 40-60% decrease in the 
umbra relative to the quiet Sun. This decrease is because 
part of the source e nerg:y is lost afte r the multiple mode 
transformations (see ICally fc Bogdanl . Il997l ). In the high- 
frequency band, however, there is an excess of the power up 
to some 40-50%. Note that the magnitude of the halo in 
our simu lations is in a good agree ment with the observed 
one (e.g. iHindman fc Brownl . [l998l ) . 

4. Discussion and Conclusions 

Our results indicate that the halo effect happens in a nat- 
ural way due to additional energy input from the high- 
frequency fast mode waves produced after their refraction. 
The halo is produced in the photospheric regions where the 
field is intermediate implying that the Alfven speed is lower 
than the sound speed. The halo is not observed at low fre- 
quencies because these waves are already reflected below 
the transformation layer. The halo is not observed in the 
umbral part of the sunspot because the refraction happens 
below the layer visible in spectral line observations. 

In our simulations, the halo is observed at the periphery 
of the strong field zone where the field is inclined by some 
30-40 degrees. For sunspots with larger field strengths than 
considered here, the level cs = va would be located deeper 
and would intersect the logrs = level further from the 
sunspot axis. Thus, for models with more intense magnetic 
field, we expect the halo will appear at larger distances 
from the sunspot axis. Also, by increasing the magnetic 
field inclination in the penumbra of the sunspot model, the 



fast to fast mode transmission will be more efficient (|Callvl . 
2006) and we can expect the magnitude of the halo to be 
somewhat larger. 

Based on our model, we can speculate regarding some 
observed properties of the halo as well as others (still un- 
detected) that may be interesting to observe in the future. 

(i) Several observation s indicate tha t the radi us of the 
halo i ncreases with height (Brown et al.L[T9 92: Bra un et al.L 
ll992l : lThomas fc Stanchfieldl . l2000D . Tliis effect can be qual- 
itatively explained by our model. As follows from e.g. FiglU 
with increasing height, the va = cs layer is located at pro- 
gressively larger distances from the sunspot axis. Observed 
in chromospheric lines, the condition necessary to detect 
the halo (i.e. that the line formation layer lies below the 
'^A = Cs layer) would be fulfilled at larger distances from 
the sunspot axis, and, in a natural way, this would produce 
halos with larger radius. We can predict that, after some 
height in the chromosphere where the whole atmosphere is 
magnetically dominated the halos should disappear. 

(ii) Our model al so explains qualitatively the observa- 
tions of iBraun fc Lin dsev (2000) who detected significant 
reflection of the high-frequency waves in active regions. 

(iii) The horizontal velocity component (not shown in 
this paper) shows a stronger magnitude of the halo effect 
as the waves propagate nearly horizontally at the heights 
where they are refracted. Thus, we suggest that the mag- 
nitude of the halo in off-center observations should be 
stronger. We are aware of only one observations of this kind 
(iToner fc Labont j . Il993l ). where apparently no change of 
the halo magnitude was detected. However, more observa- 
tions are required to confirm/discard this conclusion. 

(iv) Since the magnitude of the mode transformation 
and reflection at the va = cc layer depends on the magnetic 
field inclination (|Callvl . l2006h . we can speculate that halos, 
when detected with high resolution observations, should 
show fine structure effects in active region penumbral fil- 
aments, being more pronounced for horizontal fields. 
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